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1 Introduction 



In a recent paper Brevik, Milton and Odintsov (T| calculated the thermal energy associated 
with several types of quantum fields living in S* 1 x S^ 1 geometries with d = 2,3. In par- 
ticular the thermal energy of a quantum neutral scalar field defined on a one-dimensional 
compact space at finite temperature (S 1 x S 1 ) was found by those authors to depend 
quadratically and also linearly on the temperature in the high-temperature limit. Their 
result reads 



where (3 is the reciprocal of the thermal equilibrium temperature T. This expression is in 
agreement with the corresponding one obtained by Kutasov and Larsen |2| (in the high- 
temperature limit) but at variance with the one obtained by Klemm et al [3|]. These last 
authors justify their result by arguing that the contribution of the zero modes becomes 
significant only when the saddle-point approximation does not hold anymore. This dis- 
agreement has arisen some doubts concerning the correctness of Brevik et aVs evaluation. 
Some other authors also strongly support the necessity to get rid of the zero mode (see, for 
instance, 0), a common practice also in e.g. zeta function regularization. However, it has 
been argued by Dowker that the correct result for the thermal energy for this particular 
arrangement of quantum field and geometry contains a positive, linearly-dependent term 
in the temperature that in the high-temperature limit cancels exactly the corresponding 
negative term in (|l]), see [|] and references therein. This same result can be inferred 
from Cardy's analysis of partition functions for conformal field theories in higher dimen- 
sions ||. Besides a fundamental problem concerning the associated entropy function, this 
disagreement may have consequences for the Verlinde-Cardy relations (?], ||. 

The purpose of this brief note is to examine such questions a little bit further. We will 
see that a generalised zeta function approach to the problem shows that: (i) there is indeed 
a contribution to the thermal energy of the spatial zero mode as calculated independently 
by Dowker and Dowker and Kirsten ||1G|| , and also Cardy ||; (ii) though the thermal 
energy is independent of the scaling mass inherent to zeta function regularization, the 
free energy and the entropy associated with the spatial zero mode depend on the scaling 
mass; and finally, (iii) that there is no sense in discarding the terms that depend on the 
scaling mass, since this would lead to a violation of the fundamental thermodynamical 
relations; in fact, these need the inclusion of a dimensionful parameter to go with (3 and 
provide the unit scale in such a way that the laws of thermodynamics be obeyed by the 
thermodynamics of the zero mode. Or, to put it in other words, it does not seem to have 
a sense to choose a particular, convenient value of the scaling mass at this stage of the 
procedure, any such value being inescapably arbitrary. The full resolution of the problem 
can only be obtained after carrying out a compulsory renormalization procedure which 




a/(3 > 1, 



(1) 
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makes contact with physical experiment. 



2 Zeta function approach to the partition function 

We begin by briefly sketching a generalised zeta function approach to the partition 
function of a quantum scalar field in thermal equilibrium with a heat bath at a fixed 
temperature 

The zeta function for conformally coupled massless scalars fields living on a compact 
<ij?-dimensional manifold at finite temperature reads 





= ^EH)" S , (2) 

A 

where d 2 := V 2 + d 2 /d 2 r, \i is a scaling mass, £ = (D — 2) /4 (D — 1) is the conformal 
parameter and R is the Ricci curvature scalar. Here D = cIr + 1 where cIr denotes the 
number of dimensions of the manifold. The eigenvalues uj\ of the Euclidean d'Alembertian 
operator are given by 

where, in principle, n = 0, ±1, ±2, . . .. Here M|/a 2 denote the eigenvalues of the Laplace 
operator on the (i^-dimensional manifold. For a bosonic neutral scalar field the partition 
function can be obtained from the zeta function through 



log Z 09) = ~C' is = 0; . (4) 



We can separate the infinite temperature sector, which corresponds to n = 0, i.e., the 
zero mode associated with the compactified Euclidean temporal dimension, from the finite 
temperature sector and at the same time isolate the spatial zero mode in the following 
way. Suppose that the spatial zero mode corresponds to £ = and the remaining £ are 
shifted to £ + 1 such that £ = 0,1,2,... with an appropriate degeneracy factor taken 
into account. For convenience, however, suppose initially that the spatial zero mode is 
absent. Then taking Eq. (|3|) into (|2|), performing a Mellin transform, and making use of 
the Poisson sum formula (or Jacobi theta function identity) 



OO 1 oo 
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with r = Ant/ '/J 2 , we have 



pOO / \ it hi 

r (s) vatt f^ Jo r (s) ^ 

x W dtt s -^ 2 e- n2p2 l u - M "' a2 \ (6) 

n=l ® 

where D e is the degeneracy factor. It is not hard to see that the spatial zero mode can 
be taken into account by considering the additional term 

cz,, u 9 -^) = £ r dtts " ,2 ^ 2,At - ^ 

V v- J r(s) v % ? ^ 1 7o 

The contribution of the spatial zero mode as it stands is divergent, the origin of this 
divergence being the simultaneous consideration of two zero modes, the temporal and 
the spatial ones. Nevertheless, we can regularise it and extract its finite contribution in 
several ways. One of this ways is to apply Eq. (|5|) and discard terms that do not depend 
on /3 or contribute to the partition function with a term which is linear in (3. We obtain 



d 2 + iR\ 2 " 
/i 2 J Utt 



Cz.m. [a; 2 S 1=2(^-1 Ch(2s). (8) 




Since (r (0) ^ 0, we can expect a scaling mass dependence of some or all thermodynamical 
quantities associated with the spatial zero mode. The derivative at s = reads 

= 0;^±^)=41ogg) CB (0) + 4C'(0). (9) 

The last term is linear in (3 and can be discarded, as we did above (its contribution can 
readily be absorbed in the regularization mass /i). 

An alternative way (much more elegant mathematically) of treating the thermal con- 
tribution of the zero mode is to proceed along the lines pioneered in |T2| , that is (in a 
nutshell) to use the ordinary formulas of zeta regularization till the end, while 'lifting' the 
zero modes (e.g. preventing them from becoming zero), both the spatial and temporal 
one in this case, with the help of two small parameters (say e and r/), which are finally 
taken to zero after the whole zeta regularization process (which uses the ordinary expres- 
sions in the absence of zero modes) has been carried through. This method provides, in 
some situations, a justification of the usual principal part prescription (see an d in 
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the present case leads also to a complete cancellation (more than that, it actually avoids 
the appearence) of the divergences that we have just discarded on physical (dimensional) 
grounds. 

Making use of the integral representation for the modified Bessel function of the third 



kind 13 



poo / ft \ i^/2 

J dxx u - 1 e-*- bx = 2(j] K v (2Vab), (10) 

for the temperature-dependent part describing the non-zero modes, and adding the zero- 
mode contribution, we have 

log Z ((3) = - A p DtMt - \ log - p D e log (l - e -W) , (n) 

where we have also used the explicit form of K\ji (z) [|13|]. Notice that in this approach 
— and also in the one described in |12| — there is no temperature dependent pole as in 
1^. The factor 2ir in the denominator of the log term in this equation is fully arbitrary (it 
can obviously be absorbed into the \i parameter). It would lack any sense to set /i equal 
to a particular value (say 1, or 2tt), as \i clearly provides the (regularization) scale. 

Equation (pTT]) can be applied to several situations concerning scalar conformally cou- 
pled fields on a <iR-dimensional manifold. The first term of Eq. ([TTD is associated with the 
zero temperature vacuum energy, the second one describes the thermal corrections due to 
the spatial zero mode (when present, if not so, this term has just to be deleted), and the 
third one is the thermal contribution due to real excitations of the conformally coupled 
fields on S d ~ x . This simple and easy to interpret formula is the result of the generalised 
zeta function approach used here. 



3 The thermodynamics associated with the zero mode 

We consider, to start with, the thermodynamics associated with the zero mode. First of 
all, let us remark that the thermal energy associated with the zero mode does not depend 
on the scaling mass, in fact 

= = = (12) 

On the contrary, the free energy and the entropy do. The free energy is 

Fz.m. 09) = -- log Z ZM . {(3) = - log , (13) 
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and the entropy reads 

Sz.m. 09) = ~P 2 ^ Zz.m. 09) = - log ^ + 1. (14) 
Since the thermal energy and the free energy are related by 

(3F ((3)=jE ((3) + C, (15) 

we see that, in order to have complete compatibility between fll2"|) and (|13D , we must 
necessarily choose C = log (/x/27r). One can also verify, quite easily, that the fundamental 
relation between the three quantities above, 

S = (3(E-F), (16) 

is actually satisfied by fll2"|), (|13|) and (|14]). These results are in complete agreement with 
the ones obtained by Cardy [§], provided that we choose fi = 2n 'in inverse length units'. 
Their physical relevance, however, is dubious because, as we have already advanced and 
clearly see now, in order to maintain internally consistent thermodynamics, the need for 
a dimensionfull parameter [i is unavoidable and setting it equal to any particular value 
at this stage (that is, setting arbitrarily a mass scale) is not justified. For this reason, in 
the following we will leave the zero mode aside (however, we will discuss the effect of its 
inclusion a bit later). 



4 The high-temperature regime 

For conformally coupled scalar fields on S d ~ x , the relevant contribution from Eq. (|Tl]) 
(skipping the zero mode) reads 

oo 

log Z (/?) = - £ (£ + I)"' 2 log (l - e-^ a ) . (17) 
e=o 

And, also from Eq. the zero temperature vacuum energy is 

1 00 1 
£o = Ya E V + l)"" 2 (< + !) = (1 - d) . (18) 

Let us consider the one-dimensional compact space for which cLr = 1 and D = d = 2. In 
this case 
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On the other hand, Eq. (??) reads 



logZ(/3) = -5>g(l-e-^/ a ), (20) 

1=1 

where we have shifted the quantum number t by one unit. Since the thermal energy is 
given by 

E((3) = ~logZ(P), (21) 



one can combine (|19|), (f2~0|) and (|2l|) to obtain the total energy, in the form 



1 If I 
~ ~2^ + a^[ eW»-l' (22) 

One can now make use of an appropriate sum formula, say Euler-MacLaurin, Abel-Plana, 
or Poisson, in one of its several versions, to evaluate the sum on the rhs of ( p2l) . As an 
alternative to Euler-Maclaurin's formula employed in |E|], let us consider here the Abel- 
Plana sum formula 

£ / W = ~\f (0) + jf dx f(x)+i f dx f {lX) e ~J}^ X) ■ (23) 
A straightforward application of this expression, with / (x) = xj {eP x l a — lj, leads to 

which is the result obtained in (l| and questioned in 0. The inclusion of the zero mode 
here would have certainly led, unambiguously, to a cancellation of the (3~ l dependence. 
But this is a very specific feature of the thermal energy, which is the only thermodynamical 
quantity that does not depend on the regularization scale of the zero mode. 

Notice that ( |22|) is more appropriate for a low-temperature representation of the en- 
ergy, and upon the use of the Abel-Plana sum formula it yields the high-temperature 
result. As is well known, this is a property of this kind of sum formulas (starting with the 
Jacobi identity), which relate a series expansion valid for say low values of the relevant 
parameter with a corresponding series expansion for large values of the parameter. In 
spite of those formulas being identities (both expansions correspond of course to the same 
function), one should be careful not to mix terms of one of the expansions with terms of 
the other one. 
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It can be verified precisely that it is the third term in the Abel-Plana formula the one 
that exactly cancels the zero temperature contribution. In fact, the contribution of this 
term reads 



e 2nx _l J Q e 2nx _\ 24 



Observe also that the result given by fl24|) is closed, in the sense that the use of the Abel- 
Plana sum formula does not allow for exponential corrections, hence it holds only in the 
very high temperature regime. 



5 The low-temperature regime 

We will now construct an alternative representation for the thermal energy. Contrary to 
what has been done in the preceding section, through most of the present one we will be 
dealing with expressions valid in the high-temperature regime, and towards the end we 
will make use of the Abel-Plana formula in order to obtain the desired low-temperature 
expansion.^] 

In order to obtain, to begin with, a high-temperature representation of the free energy, 
we make use of the Poisson sum formula, in the form 

V) / {n) = - J -±-L + / di/(i) + 2V dxf(x) cos (2imx) , (26) 

n=l 2 J ° n=l J ° 



where here / (x) = log (l — e P x l a \ . A straightforward application of (|26|) to (|20| ) yields 

f (0) 00 
logZ(/?) = ^-J(0)-2£/(n), (27) 

/ 71=1 

where 

I(n):= dx f (x) cos (2imx) , n = 0, 1, 2, 3, . . . . (28) 
Jo 

The integrals defined by I (n) can be readily evaluated with the help of formulas to be 
found in |I3| , and the result is 

, , ^ f(0) ir 2 a f3 1 ~ 1 , f2iT 2 na\ . , 



1 Both the high and low temperature ex pan sions can be derived, in a completely equivalent but more 
direct way, from the representation of Eq. (|2^ ) in terms of the normalized Eisenstein series E2 (r) [§, . 
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The contribution due to / (0) can be explicitly calculated as follows. Expanding the log 
as usual we can write 

f(x)= daI»C(a + l) ^ (30) 

2m Jc-ioc \ a J 

where we have made use of the representation 

/■c+ioo 

/ daT{a)s~ a 1 s, c real, \s\ < 1, c > 1. (31) 

J c—ioo 



] rc+ioo 

2id 



In the limit x — > the finite contribution of / (x) is determined by the double pole of the 
integrand in Eq. ( j30| ) at a = 0, 

lim /W = /(0) = log(Aj. (32) 
On the other hand, we can also make use of the expansion 



coth (x) = 1 + 2 e~ 2px , x>0. (33) 
P =i 

As a consequence, after discarding an unphysical (formally divergent) constant we end up 
with 

^^■J^^-^fW. (34) 
It follows that the thermal part of the free energy in this representation is given by 

P{p) = -iiogzoa) 

6^ 2/? l0 H2vraJ + 24a /? ^ n 6 ' (35) 

As in the case of Sect. 2, an alternative calculation can be performed using standard 
zeta function methods ab initio, namely, the transform (|3lD already in the expansion 
of the starting expression (??). The final result obtained is exactly the same, but the 
calculation is very much shortened and (as before) it avoids at every step the formation 
of divergences.^ 

2 We owe K. Kirsten this observation and the correction of an error in our previous result. 
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Notice that our result, Eq. (|35|), has a log term and therefore does not fulfill the 
general relations given in @, namely: that the high temperature expansion of the free 
energy must have the general form 



2ira \ d ( 2-ira \ d 2 1 2ira * 



-F{(3)a = aA — \ + a d - 1 aA — \ + "-- + a o( — J + O (V 47r a/p ) , (36) 
the coefficients of this expansion satisfying the relation 

d/2 

£(-l) fc (2fc-l)a 2fe = 0. (37) 

fc=0 

If not for the log term, these relations would be fulfiled. Observe, however, that there 
would not be an immediate cancellation of the log term if we had included the zero mode 
contribution. In fact, the matching of the zero mode term, log with the one, with 

opposite sign, coming from the high-T limit, -log (Jj^) , is not a triviality. On the contrary, 
it has a very deep meaning: it takes place if and only if we prescribe that the regulating 
mass-dimensional parameter be set exactly equal to the inverse of the compactification 
length, fi — 1/a. This may certainly be viewed as a most natural choice and a physically 
meaningful one — being 1/a the only caracteristic dimension in our model — but even 
then, this is no substitute for the pertinent renormalization procedure. 
The total free energy is 

F (fl = E + F W = -g - ± log (JL) - 1 £ I e — /, (38) 
The total energy reads 

E (fl = | (PF W) = g - ± - 1 1 V e-—<°, (39) 

and can be summed over n to yield the high temperature representation 

n 2 a 1 47r 2 a ^ p 



If we now, as before, make use of the Abel-Plana sum formula (|23|) , we obtain 

E ^) = ~^ (41) 

Notice that, as anounced above, we have here started from a high-temperature represen- 
tation for the energy and now, upon the use of the Abel-Plana sum formula, we obtain the 
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energy in the zero-temperature limit, a procedure that can be considered the reciprocal of 
the one employed in Ref. Notice also that in the Abel-Plana formula, the first term 
leads to the zero temperature term, the second cancels out the linear contribution and 
the third one cancels out the Stefan-Boltzmann term. 



6 Final remarks 

In this brief analysis, we have focused on the possible contribution of the spatial zero 
mode to the thermal energy of a scalar field in a S 1 x S* 1 geometry. We have shown that, 
though the thermal energy does not depend on the scaling mass //, the partition function, 
the free energy and the entropy do depend on the choice of value for this parameter. The 
important obstruction to be circumvented is to satisfy the third law of thermodynamics 
for, as it stands, the entropy (|H|) associated with the spatial zero mode diverges in the 



zero temperature limit, a fact also noticed in 0. It was argued in || that the scaling 
mass can be redefined. Unfortunately, any tentative of redefining the scaling mass, that 
is, of eliminating the scaling mass dependence from (PH), will lead to clear contradictions 
among the thermodynamical relations linking the basic quantities (|T^), (|13|) and (fl-il) . In 



fact, the scaling mass parameter cannot be made disappear. At most, it can be traded 
for a characteristic mass or inverse length of the system considered (here 1/a, the inverse 
compactification radius). A subsequent renormalization procedure seems inescapable, 
which makes contact in the end with physical reality. 

Let us point out that, although we have here considered a very simple model, in order 
to keep all terms under easy control and render the argument clear, the results obtained 
are valid under much more general circumstances, affecting very general physical systems 
that develop zero modes. The treatment has been that of (a natural extension of) the 
zeta function regularization prescription. 

To sum up, only if one restricts the analysis to the thermal energy — and nothing 
else — can one avoid the situation here presented, since this is the only thermodynamical 
quantity that does not depend on the scaling mass. It would seem, thus, that coherence 
forces one to set the zero mode and its thermodynamics aside. Notice also that the thermal 
energy, the entropy and the free energy of the zero mode do not explicitly depend on the 
radius (or characteristic compactification length) of the manifold. This means that (at 
least at the one-loop level) their contribution to any isothermal process — e.g., the only 
ones permitted if the system is in thermal equilibrium with a heath bath — will not be 
physically observable. But this consideration would be completely reversed if the necessity 
for setting in fact fM—l/a would prevail in the end (as one of the factions in the already 
mentioned controversy tend to support), and this would open the door to an experimental 
manifestation of the zero mode and to the possibility of answering in the affirmative the 
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important question of the controverted physical relevance of these zero modes. In that 
sense, a fine experimental validation of Eq. (|24]) would be crucial (presence or absence of 
the linear dependence). For what we have seen, this would imply, necessarily, the breaking 
of the thermodynamical relations and the necessity of reformulating the third principle in 
the quantum field theory domain. 
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